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1.  INTRODUCTION 

One  of  the  major  issues  in  the  hindcasting  and 
forecasting  of  winds  and  waves  is  the  difficulty  of 
validation  and  verification.  While  comparisons  wi* 
point  measurements  from  discrete  and  sparsely 
distributed  wave  buoys  provide  some  measure  of 
statistical  confidence,  the  spatial  distribution  of  the 
modeled  wind  and  wave  fields  cannot  be  easily 

assessed.  . 

Remote  sensing  from  space  provides  a  synoptic 

view  of  the  ocean  wind  and  wave  fields.  For  example 
wind  speed  and  significant  wave  height  are  standard 
output  of  spaceborae  altimeters  such  as 

TOPEX/POSEIDON  (hereafter  referred  to  as  lUl'tA 
for  brevity).  Comparisons  of  the  altimeter  measured 
wind  speed  and  wave  height  with  surface  buoy  data 
have  shown  very  positive  agreement  (e.g.,  Ebuchi  and 
Kawamura  1994;  Freilich  and  Challenor  1994,  Gower 
1996;  Hwang  et  al.  1997).  With  an  along-track 
resolution  of  7  km,  the  spatial  resolution  of  the 
spacebome  altimeter  is  comparable  with  that  of  the 
numerical  models  used  for  regional  simulations.  In  the 
following,  we  present  the  results  of  a  companson  study 
of  WAM  wave  modeling  of  the  Yellow  and  East  China 
Seas  (YES)  with  two  of  the  TOPEX  tracks  m  the 
region.  In  the  next  section  the  accuracy  of  TOPEX 
altimeter  wind  and  wave  measurements  is  summarized 
based  on  earlier  studies  of  TOPEX  and  surface  buoy 
comparisons.  Section  3  describes  the  YES  data  sets  and 
the  background  information  about  the  numerical 
modeling  and  satellite  tracks  in  the  comparison  region. 
Section  4  presents  the  results  of  the  companson  in 
terms  of  statistics  such  as  bias  (5),  rms  difference  (A), 
regression  coefficient  (c  and  Cy,  to  be  further  discussed 
later)  and  correlation  coefficient  (R).  Since  altimeter 
remote  sensing  provides  a  spatial  coverage  of  the  win 
and  wave  fields  along  transects,  we  will  explore  the  use 
of  such  information  for  the  validation  of  wave  he^t 
distribution  from  a  numerical  wave  model  output.  e 
conclusions  and  summary  are  presented  in  the  last 
section  of  the  paper. 

2.  ALTIMETER  MEASUREMENT  OF  WIND  AND 
WAVES 

Extensive  studies  of  altimeter  measurements  of 
wind  and  wave  parameters  have  been  publishqj  (e  g-> 


Brown  et  al.  1981;  Chelton  and  WenU  1986;  Wu  1992; 
Witter  and  Chelton  1994;  Freilich  and  Challenor  1994; 
Ebuchi  and  Kawamura  1994;  Gower  1996;  Hwang  et 
al  1997  and  the  references  therein).  These  studies 
have  shown  consistently  that  the  significant  wave 
height  derived  from  the  altimeter  is  essentially 
identical  to  that  measured  by  ocean  surface  buoys. 
Particularly,  if  the  distances  between  the  buoy  location 
and  the  satellite  footprints  are  within  10  km, 
difference  of  the  two  sets  of  data  is  approximately  0. 
m.  One  exception  is  in  the  coastal  region,  where  large 
disagreement  often  occurs.  The  cause  of  the 
discrepancy  may  be  due  to  the  land  effect  on  radar 
backscatter  or  the  greater  spatial  inhomogene^  of  the 
coastal  wave  field  within  the  radar  footpnnt.  There  are 
no  established  rules  relating  the  detenoration  of 
wind/wave  data  to  the  land  proximity.  Factors  effecting 
the  distance  of  data  deterioration  include  the  satellite 
approach  (from  land  or  from  sea),  the  local  bathymetry 
gradient,  and  the  plane  geometry  of  the  subject  location 
and  shoaling  factor.  From  exaimnmg  the  satellite 
output,  it  is  estimated  that  if  the  satellite  approach  is 
from  land  to  sea,  three  to  four  footprint  lengths  (20  to 
30  km)  are  needed  for  the  signal  to  adjust  to  the  change 
of  backscattering  characteristics  between  land  and 
water.  For  the  sea  to  land  approach,  this  distance  is 
reduced  to  one  footprint  (7  km),  but  odier  factors 
mentioned  above  may  dictate  the  application  distance 
from  land.  Since  the  issue  concemmg  land  proxmuty 
has  not  been  thoroughly  investigated,  we  will  exclude 
data  within  50  km  of  land  in  the  following  companson. 

Derivation  of  wind  speed  from  microwave  radars  is 
based  on  the  correlation  of  surface  roughness  and  radar 
backscatter  intensity  (the  normalized  radar  backscatter 
cross  section,  Ob)-  the  spacebome  altimeters, 

specular  reflection  is  the  primary  mechanism  of  radar 
backscatter.  In  this  mode,  roughness  on  the  surface 
causes  incidentradar  waves  to  diffuse  and  scatter  away 
from  the  aperture  of  radar  reception,  and  the  smoother 
the  surface  the  higher  the  backscatter  is  expected.  On 
the  ocean  surface,  short  ocean  waves  are  the  primary 
contributor  to  the  surface  roughness.  Because  these 
short  water  waves  are  wind  generated,  higher  Oo 
corresponds  to  lower  wind  speed  and  Y"", 
the  past  two  decades,  the  operational  algorithms  to 
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retrieve  wind  speed  from  the  altimeter  backscatter 
cross  section,  however,  have  been  based  on  empirical 
formulae  established  On  collocated  databases  of 
altimeter  and  surface  measurements  (e.g..  Brown  et  al. 
1981;  Chelton  and  Wentz  1986;  Witter  and  Chelton 
1991)  or  on  the  statistical  distributions  of  wind  speeds 
and  Oo  separately  (Freilich  and  Challenor  1994).  These 
comparisons  indicate  that  despite  the  continuous 
improvement  of  altimeter  hardware  and  software,  the 
progress  of  wind  speed  retrieval  from  various 
algorithms  developed  over  the  last  20  years  only  shows 
marginal  improvement.  For  example,  Freilich  and 
Qiallenor  (1994)  report  a  very  comprehensive  analysis 
of  eight  algorithms.  Results  in  terms  of  the  mean  error 
(bias),  root  mean  square  eiror,  standard  deviation,  wind 
speed  error  trend  and  pseudo-wave  age  error  trend  are 
tabulated.  The  differences  in  all  eight  algorithms  are 
relatively  minor.  In  terms  of  mean  square  errors  and 
standard  deviations,  7  of  the  8  algorithms  produced 
less  than  10%  difference  in  these  two  error  parameters 
(1.60  to  1.75  m/s  rms  error,  1.58  to  1.72  m/s  standard 
deviation).  When  applied  to  a  regional  scale,  the  error 
statistics  are  slightly  better.  For  example,  comparison 
of  3.2  years  of  TOPEX  data  with  4  buoys  in  the  Gulf  of 
Mexico  suggests  that  when  the  separation  distance 
beriveen  the  TOPEX  footprint  and  the  buoy  location  is 
less  than  10  km,  the  rms  difference  is  approximately  1.2 
m/s  (Hwang  et  al.  1997).  Figure  1  shows  an  example  of 
the  comparison  of  significant  wave  height  and  wind 
speed  measured  by  the  TOPEX  Ku-band  altimeter  and 
the  surface  buoys  in  the  Gulf  of  Mexico. 


Figure  1.  Comparison  of  the  significant  wave  height  and 
wind  speed  measured  by  the  TOPEX  altimeter  and  the 
ocean  buoys  in  the  Gulf  of  Mexico  (Hwang  et  al.  1997). 


The  regression  coefficient,  rms  difference  and 


correlation  coefficient  for  the  wave  height  are  0.98, 
0.13  m  and  0.97,  respectively;  the  corresponding 
statistics  for  the  wind  speed  are  1.0,  1.2  m/s  and  0.85 
for  this  case. 

3.  YELLOW  AND  EAST  CHINA  SEAS  DATA  SET 

As  illustrated  in  the  last  section,  the  altimeter 
output  of  winds  and  waves  represents  a  good  source  of 
wind/wave  measurements  for  the  study  of  regional 
wind/wave  climate.  The  data  can  also  be  used  for 
verification  and  validation  of  wave  model  output  and 
offer  insight  into  the  spatial  distribution,  which  can  not 
be  verified  easily  with  point  measurements.  This  is 
especially  valuable  in  regions  where  in-situ  data  are 
either  nonexistent  or  inaccessible. 

The  WAM  model  (Cycle  4)  (WAMDl  Group  1988; 
Janssen  1991)  was  recently  applied  to  the  Yellow  and 
East  China  Seas  to  hindcast  the  wave  climatology  for 
1994.  The  wave  hindcast  setup  consisted  of  a  global  1“ 
Lat/Lon  grid  to  generate  input  wave  boundary 
conditions  and  a  regional  0.25  °  Lat/Lon  grid  including 
the  Yellow  Sea,  East  China  Sea,  and  the^Sea  of  Japan 
shown  in  Figure  2. 


Figure  2.  The  model  domain  of  the  WAM  simulation  for 
the  year  1994  used  in  this  simulation. 


A  global  wind  product  (Wobus  and  Kalnay  1995) 
from  the  National  Centers  for  Environmental 
Prediction  (NCEP,  formerly  the  National 
Meteorological  Center)  was  used  to  drive  the  global 
and  regional  wave  model  grids.  The  NCEP  winds, 
which,  are  provided  at  approximately  0.94°  Lat/Lon 
resolution,  were  bilinearly  interpolated  to  a  1° 
resolution  for  the  global  hindcast  and  a  0.25° 
resolution  for  the  regional  hindcast.  The  global 
hindcast  was  run  with  a  temporal  resolution  of  6  hours, 
which  is  the  resolution  provided  by  NCEP.  The  NCEP 
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winds  were  linearly  interpolated  to  a  3  hour  resolution 
for  the  regional  hindcast. 

Several  researchers  have  suggested  improvements 
to  WAM  (Tolman  1992;  Lin  and  Huang  1996;  Bender 
1996).  Bratos  (1997)  compared  deepwater  WAM 
results  to  National  Data  Buoy  Center  (NDBC)  wave 
measurements  along  the  US  Atlantic  coast  using  high 
quality  winds  (Cardone  1992)  for  five  extratropical 
storms  which  included  a  variety  of  conditions  ranging 
from  extreme  events  to  more  moderate  and  variable 
events.  One  persistent  tendency  found'  in  the  WAM 
results  was  a  low  bias  significant  wave  height  (SWH), 
Figure  3  shows  the  relationship  between  wind 
speed  bias  and  SWH  bias  for  all  buoy  locations  and 
storm  events  considered.  For  a  well  behaved  model  the 
SWH  bias  should  correspond  to  the  wind  speed  bias. 
The  plot  shows  a  tendency  for  WAM  SWH  to  be 
biased  low  (0.0  m  to  -0.5  m)  when  the  wind  speed  is 
biased  high  (0.0  m/s  to  1.75  m/s). 

One  of  the  standard  WAM  output  files  includes  the 
SWH,  friction  velocity  w*  ,  and  drag  coefficient  at 

every  grid  point.  The  model  wind  speed  at  a  10  m 
elevation  is  given  by  .  The  SWH  and 

w'ind  speed  model  parameters  were  extracted  along  the 
two  TOPEX  tracks  for  comparison.  The  first  track  is 
number  69,  a  descending  track  that  runs  along  the 
central  axis  of  the  YES  region  (Figure  4),  and  the 
second  track  is  number  26,  which  is  an  ascending  track 
cuning  across  the  entrance  region. 
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Fig.  3.  Example  of  WAM  output  compared  with  buoy 
measurements  along  the  U.  S.  Atlantic  coast  [Bratos 
1997]. 

The  extraction  of  the  model  output  consisted  of 
interpolation  of  the  model  parameters  to  the  satellite 
track  coordinates  using  two  subroutines  within  the 
GMT-System  (General  Mapping  Tools)  (Wessel  and 
Smith  1991).  The  subroutine  surface  (  Smith  and 
Wessel  1990)  creates  a  binary  grid  file  using  an 


adjustable  tension  continuous  curvature  algorithm.  The 
subroutine  grdtrack  samples  this  2-D  binary  grid  file 
along  tracklines  26  and  69  using  bicubic  inteipolation. 
The  numerical  output  and  altimeter  data  are  then 
merged  together  by  linear  interpolation  using  the 
latitude  coordinate  to  a  uniform  spacing  of  0.05  degree 
in  latitude.  The  spacing  between  neighboring  data 
points  in  the  merged  data  set  is  between  6.0  to  6.3  km 
along  track.  The  distances  between  the  model  data 
positions  and  the  interpolated  altimeter  footprints  are 
between  0.4  to  0.6  km.  Over  the  one-year  (1994)  model 
run,  32  cycles  are  extracted  for  each  track  (4  of  the 
missing  cycles  occur  at  Julian  days  63,  172,  311  and 
361  for  track  26,  and  76,  176,  315  and  324  for  track 
69.).  The  time  differences  between  the  model  output 
and  TOPEX  measurement  are  distributed 
approximately  evenly  between  0  to  1.5  hours. 


Longitudi 

Figure  4.  A  map  of  the  Yellow  and  East  China  Seas  (YES) 
region  for  this  comparison  study.  TOPEX  tracks  are 
shown  in  dotted  curves.  Tracks  69  (from  NW  to  SE  going 
through  the  axis  of  the  region)  and  26  (from  SW  to  NE 
across  the  YES  entrance)  are  used  for  comparing  with  the 
WAM  model  simulation  for  the  year  1994. 

4.  COMPARISON  AND  DISCUSSIONS 
4.L  Average  Comparison 

We  fust  study  the  wave  and  wind  climate  derived 
from  TOPEX  measurements  and  WAM  output.  The 
annual  averages  Slong  the  two  tracks  are  listed  in  Table 
1(A).  The  statistical  distribution  of  the  winds  and 
waves  in  four  different  regions  designated  as  are 
shown  in  Figure  5  (top  half),  in  which  tbe  histogram  of 
wind  speeds  and  wave  heights  are  displayed.  These 
four  regions  are  designated  as  D9,  S3,  M3  and  N3  for  9 
degrees  (latitude),  Southern  3  degrees.  Middle  3 
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^grees  and  Northern  3  degrees  along  the  two  tracks 
limiting  latitudes  are  tabulated  in  the  first  column 

that  tt,  ^  tabulated  statistics,  it  is  found 

Aat  the  average  properties  of  wind  speed  and  wave 
height  between  the  numerical  simulation  and  the 
remote  sensmg  measurement  are  in  excellent 
agreement.  Both  the  mean  and  standard  deviation  of 
ese  average  quantities  are  mostly  in  agreement  within 
10  percent.  Both  TOPEX  and  WAM  show  that  the 
average  wave  height  decreases  persistently  northward 
Along  the  central  axis  of  YES  (track  69)  the  3-degree 
average  wave  height  is  reduced  from  1.70  m  at  29®  to 

m  at  32  and  1.16  m  at  35®  based  on  TOPEX 
measurements. 


Figure  5.  Comparison  of  the  histograms  of  wave  heights 
and  wind  speeds  from  TOPEX  measurements  and  WAM 
fWAM^r/  arranged  as  Hg  (TOPEX),  Hj 

‘0  ‘be 

urth  column.  The  upper  four  rows  of  the  plots  are 
N  3®  (S3),  Middle  3®  (M3), 

Inw  f°  ^  *^®*P®^‘’''ely!  along  track  69.  The 

lower  four  rows  are  the  histograms  of  the  wind  and  wave 
parameters  m  the  9®  latitude  along  track  69,  and  over  the 
first  to  the  fourth  90-day  period  of  1994  data  sets. 

The  results  show  that  the  WAM  hindcast  of  the 
wave  heights  are  slightly  higher,  by  9®/o,  3®/o  and  4®/o 

(S3,  M3  and  N3).  The  average  wave  height  along  track 
26  shows  a  sumlar  trend,  with  1.82  m  at  27®,  1.63  m  at 
30  and  1.51  m  at  33®.  The  WAM  output  again 
oveipredicts  the  wave  height  by  9%,  11%  and  2®/o, 


respectively.  TTie  distribution  of  the  average  wave 
heights  can  be  measured  in  terms  of  the  standard 
deviation  (Table  1).  In  the  along-axis  direction  (track 
69),  the  dimensionless  standard  deviations  (normalized 
by  the  mean)  are  0.45,  0.65  and  0.58  from  South  to 
No^  based  on  the  TOPEX  data,  and  0.37,  0.55  and 
0.74  respectively,  based  on  the  WAM  output.  The 
Standard  deviations  along  track  26  of  TOPEX  and 

(fO-54.  0.47.  0.69] 

xiiaIS  T  “  to  [0.52,  0.53.  0.61]  of 

average  predicted  by 
WAM  model  agrees  very  well  with  the  TOPEX 
measurement.  One  of  the  major  sources  of 
discrepancies  is  the  wind  speed  or  wind  stress  used  to 
dnve  the  model.  As  shown  in  Table  1.  the  average 
wmd  speed  used  in  WAM  is  in  general  higher  than  the 

S3),  the  WAM  average  wmd  is  4  to  19®/o  higher  than 
the  corresponding  TOPEX  data. 

The  comparison  statisrics  of  the  seasonal  average 

between  TOPEX  and 
WAM  are  also  listed  in  Table  1(B).  The  histograms  of 
the  wave  height  and  wind  speed  distributions  are 
shovn  m  Figure  5  (lower  half).  These  statistics  are 
based  on  *e  data  along  the  frill  9-degree  tracks 
discussed  above.  The  statistics  over  a  smaller  or  larger 
distance  can  be  generated  also.  Two  of  the  notable 
features  of  the  seasonal  average  are  (1)  the  lower 
average  occurs  in  the  second  quarter,  and  (2)  the 
largest  variation  as  represented  by  the  standard 
deviation  occurs  in  the  third  quarter  (Table  IB)  These 
feanires  are  consistent  with  the  monsoon  climate  of  the 
region,  with  predominantly  northwesterly  to  northerly 
winds  m  the  wmter,  and  southwesterly  to  southeasterly 
wmds  m  the  summer.  In  the  transition  seasons,  the 
wmds  fluctuate  between  the  two  dominant  weather 
systems  (e.g.,  Wang  and  Aubery  1987).  The  wind  and 
wave  properties  are  further  complicated  by  the 
typhoons,  which  generally  occur  from  July  to  October, 
and  explam  the  large  variation  (standard  deviation)  of 
the  seasonal  statistics  di^g  the  third  quarter.  The 
agreement  of  the  mean  properties  derived  from  WAM 
ou^ut  and  TOPEX  measurement  is  generally  good,  the 

nT  ,  I  (WAMH-OPEX)  ranges  between 

0.94  to  1.19.  The  agreement  in  the  standard  deviation 
^tistics  IS  also  similar  to  that  of  the  seasonal  means. 

Ihe  larger  difference  in  the  seasonal  means  (up  to  19®/o 
higher  as  compared  to  ll»/o  maximum  in  the  annual 
means)  is  closely  related  to  the  wind  condition  used  in  ■ 
dnvmg  the  wave  model.  The  maximum  difference  of 

‘S  36®/o 

1.36;  Track  26,  Q2)  as  compared  to  19%  in 
the  annual  means  (7.76/6.50=1.19;  Track  26.  M3). 
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We  then  compare  the  statistical  correlation  of  the 
WAM  and  TOPEX  along  each  track.  The  region  of 
comparison,  defmed  by  the  latitude  distance  along  the 
tracks,  is  varied  from  less  than  1  degree  to  9  degrees. 
The  mean  properties  averaged  within  the  specified 
region  will  be  presented  below.  The  statistic 
parameters  of  the  linear  regression  coefficients  Cy  and 
c,  the  bias,  the  rms  difference  and  the  correlation 
coefficient  are  computed  and  summarized  in  Table  2. 
The  coefficient  c  is  calculated  by  minimizing  the 
orthogonal  distances  between  data  points  to  the  linear 
regression  curve.  In  the  scatter  plot  of  WAM  data  vs. 
TOPEX  data  with  the  intercept  forced  to  zero,  the 
regression  fitting  can  be  expressed  asy  =  cx,  where  y  is 
the  WAM  wave  height  (wind  speed)  and  x  is  the 
TOPEX  wave  height  (wind  speed).  The  formula  of  the 
orthogonal  regression  coefficient  is  described  in  Bauer 
et  al.  (1992)  and  also  in  Romeiser  (1993), 


c  =  tanj 


1  -1 
-tan 
2 


<2;cy) 


(1) 

where  angle  brackets  denote  ensemble  average.  A 
slight  modification  to  consider  the  condition  that  the 

magnitude  of  ^y^^  exceeds  is  given  by  Hwang 
et  al.  (1997) 


tan“ 


<2xy) 


+  n 


(lb) 

The  regression  coefficient  based  on  minimizing  the 
orthogonal  distance  differs  from  those  obtained  from 
minimizing  the  mean  square  distances  in  the  y- 
direction,  Cy,  or  in  the  x-direction,  Cx’. 

c  =M  ,  (2) 

'  (*■) 
and 


When  the  correlation  of  the  two  data  sets  are  high, 
these  three  coefficients  become  indistinguishable  (see 
discussions  in  Bauer  et  al.  1992). 

The  average  wave  height  and  wind  speed  of  the 
WAM  output  and  TOPEX  data  are  in  good  agreement 
from  cycle  to  cycle.  Figure  6  shows  an  example  of  the 
results.  For  the  wave  height,  the  rms  difference  for  the 
9-degree  average  is  approximately  0.3  m,  with  a 
correlation  coefficient  of  0.93  for  both  tracks  (Table  1). 
The  rms  difference  increases  to  approximately  0.4  m 


when  the  averaging  distance  decreases  and  the 
correlation  coefficient  also  decreases,  possibly  due  to  a 
decrease  of  data  points  in  the  averaging  process.  For 
the  wind  speed,  the  rms  difference  is  1 .92  m/s  for  track 
26  and  1.33  m/s  for  track  69,  with  the  correlation 
coefficients  of  0.83  and  0.91,  respectively  for  the  9- 
degree  average.  The  statistics  for  the  3-degree 
(Southern,  Middle,  and  Northern)  averages  are  also 
tabulated  in  Table  2.  In  general,  the  agreement  in  wave 
height  between  TOPEX  and  WAM  are  good, 
considering  that  the  wind  input  used  for  driving  the 
model  is  considerable  difference  from  the  TOPEX 
measurements. 


Figure  6.  An  example  of  the  mean  wave  height  (a)  and 
wind  speed  (b)  averaged  along  the  track  of  various 
distances.  The  example  shown  is  the  3-degree  average 
(M3,  from  30  to  33  degrees  latitude).  O:  TOPEX,  x  = 
WAM.  The  statistics  of  comparison  are  listed  in  Table  1. 

4,  2.  Spatial  Distribution 

One  of  the  most  significant  advantages  of  remote 
sensing  measurements  is  the  ability  to  provide  a  quasi- 
instantaneous  meastirement  of  the  spatial  distribution 
of  the  wave  heights  and  wind  speeds.  For  example,  the 
9®  latitude  transects  of  tracks  26  and  69  were 
completed  in  168  and  190  s,  respectively.  The  high 
spatial  density  of  the  measurements  represents  a  unique 
source  of  wind  and  wave  data  that  can  be  used  for 
model  comparison. 

Figure  7  illustrates  the  cycle-by-cycle  comparison 
of  wave  heights  (top  half)  and  wind  speeds  (bottom 
half)  derived  from  WAM  and  TOPEX  track  69.  For 
this  along-axis  track,  the  WAM  wave  height  output 
tends  to  be  biased  high  in  the  first  six  to  nine  months  of 
the  year  and  the  bias  trend  reverses  toward  the  winter 
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scatter  is  quite  large  (Figure 

the  '  with 

month,"^  “'=gative  bias  appearing  in  the  first  four 
months,  and  positive  the  rest  of  the  year.  While  there 
are  apparent  sMarities  of  wave  bias  and  w^d  Was 
the  scattemg  of  data  points  is  large.  The  average 
normalized  bias  relative  to  the  mean  wave  height  or 


sa 


^rto?hIin comparison  of  the  wave  height 

(  P  a'O  and  wind  speed  (bottom  halo  from  TOPEX 

(dots)  and  WAM  (solid  curve). 

normalized  (A(//^y<//^ 
mQy<Uio>)  differences  of  wave  heights  and  wWd 
pee*  ploned  ip  Figure  ge  „a  g,  foryeckfld  e"d 
69  respectively.  He  nm  differences  relative  to  the 
mean  wave  height  are  mostly  less  than  0  5  with  1 
everage  value  of  0.27  for  nach  26  and  0.24  VrLh 

0  55  and  0  77  /'"m'"'  *'  '"rrespondmg  numbers  ate 

ana  0.77  for  the  two  tracks. 

and  coefficient,  c.  for  tracks  26 

d  69  is  shown  in  Figures  8e  and  8f  As  can  be  seen 
from  the  plots,  there  are  only  small  members  of  cases 


wiA  0  9  ^  c  ^  1.1.  Tbe  WAM  output  tends  to  over 
estimate  the  wave  height,  especially  in  the  sunder 
months.  Many  of  the  wave  height  over  estimationTTan 

userasT  °''u  Ae  wind  speed 

sed  as  the  input  to  the  wave  model.  Exceptions  are 

found  m  veiy  low  wind  conditions,  where  the  wave 
heights  are  mamly  dominated  by  the  swell  components 
that  are  unrelated  to  the  local  wind  condition. 


Of  comparing  spatial  distributions  of 

WAM  output  with  TOPEX  measurements.  Bias 

normalized  by  the  maximum  for  track  26  (a)  and  69  (b)- 
bias  normalized  by  the  mean  wave  height  or  meL  wLd 
speed  for  track  26  (c)  and  69  (d);  rms  difference 
normalized  by  the  maximum  for  track  26  (e)  and  69  (fl- 
rms  dif^rence  normalized  by  the  mean  w  ve  height  or’ 
mean  wind  speed  for  track  26  (g)  and  69  (h);  regres  i^n 
coefficent  for  track  26  (i)  and  69  G);  and  corrflation 
coefficient  for  track  26  (k)  and  69  (I).  o;  wave  height  x- 
wind  speed.  (M3  averages). 

TOPEX  data  show  the  largest  variation  (Figure  8g  and 
8h)  compared  to  the  other  statistical  parameters^  For 
the  a  ong-axis  track  69.  7  out  of  32  cycles  show 
negative  conelation  coefficient,  indica  Jg  L,  Z 
opposite  trends  of  wave  height  distribution  along  21 
altimeter  track  obtamed  from  WAM  modelinj  as 

SieTeolH‘°  In  addition  to 

the  negative  correlation  cases,  9  more  cycles  have 

onrhaTfir^^'"’^  so  altogether 

one  half  of  the  cases  examined  for  track  69  show  poor 

correlation  of  WAM  output  with  alJete 

s^rerr-tb 

statistics,  for  the  wmd  speed,  there  are  7  cycles  of 
negative  correlation  and  9  additional  cases  of  low 
correlation  coefficients.  While  many  cases  of  poor 

coiJ!  !^"  oonespond  to  those  of  poor 
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Because  the  correlation  coefficient  is  primarily  a 
measure  of  proportionality  of  two  variables  (for 
example,  wheny  =  ex'*,  where  c  and  n  arc  constant,  the 
correlation  coefficient  always  equals  to  unity)  the 
correlation  coefficient  alone  does  not  indicate  the 
quantitative  agreement  of  the  two  variables.  The 
evaluation  of  data  agreement  needs  to  take  into  accoimt 
other  statistics  such  as  the  bias,  nns  difference  and/or 
regression  coefficient. 

The  cycle-by-cycle  statistics  of  track  26  are  not 
very  different  from  those  of  track  69.  For  the  wave 
height,  there  are  4  cases  of  negative  correlation 
coefficient  and  8  cases  with  0  <  /?  <0.75.  For  the  wind 
speed,  there  are  9  negative  correlation  and  10  cases 
with  Q<R  <0.75.  Taking  the  statistics  of  the  two  tracks 
together,  only  in  about  one  half  of  the  cases  are  WAM 
output  and  TOPEX  measurements  in  good  agreement. 
The  above  conclusions  are  based  on  analyzing  32 
cycles  of  2  tracks  over  a  3-degree  latitude  stretch  of  the 
satellite  tracks.  Taking  a  different  stretch,  the  statistics 
are  not  significantly  different  from  the  above  results. 

5.  SUMMARY  AND  CONCLUSIONS 

The  wind  speeds  and  wave  heights  measured  by 
satellite  altimeters  have  been  compared  with  ocean 
buoy  measurements  with  excellent  agreement.  They 
represent  a  good  source  of  environmental  data  for  the 
study  of  global  and  regional  wave  conditions.  In  this 
paper,  we  compare  TOPEX  altimeter  measurements 
with  WAM  model  output  in  the  region  of  Yellow  and 
East  China  Seas.  Two  of  the  satellite  tracks  passing 
through  the  region  are  selected  for  comparison.  Track 
69  passes  through  the  central  axis  of  YES  and  track  26 
passes  across  the  entrance  to  the  enclosed  seas.  The 
comparisons  are  processed  into  two  categories:  (1) 
comparing  the  time  series  for  a  given  location,  and  (2) 
comparing  with  the  spatial  distribution  for  a  given 
time. 

The  first  category  is  further  divided  into  (a)  average 
over  a  region  from  3  to  9  degrees  latitude  along  a  track, 
and  (b)  seasonal  average  along  a  track.  These  average 
quantities  from  TOPEX  and  WAM  are  in  good 
agreement  in  terms  of  mean  and  standard  deviation  of 
the  statistics  (Table  1).  The  difference  between 
altimeter  measurements  and  numerical  simulations  is 
approximately  10%  with  WAM  overestimating  most  of 
the  time.  The  median  value  of  the  rms  differences  is 
0.4  m  with  a  correlation  coefficient  of  -0.9  (Table  2). 
These  figures  are  compared  to  0.15  m  and  0.97, 
respectively  from  TOPEX/buoy  comparisons 
performed  in  the  Gulf  of  Mexico  region  for 
comparable  spatial  and  temporal  separations  (Hwang  et 
al.  1997).  The  rms  difference  in  wind  speed  in  the 
YES/WAM  comparison  is  -2  m/s  with  a  correlation 


coefficient  of  *-0.83.  The  corresponding  values  from 
TOPEX/buoy  comparisons  in  the  Gulf  of  Mexico 
region  are  1.2  m/s  and  0.9,  respectively. 

The  comparison  of  spatial  distributions  of  WAM 
simulations  and  TOPEX  measurements  is  performed 
with  32  cycles  each  from  the  two  tracks  over  3  to  9 
degrees  latitude  (Figure  7  shows  the  results  of  9-degTee 
processing).  The  quality  of  agreement  between 
numerical  simulations  and  altimeter  measurements 
varies  significantly  from  cycle  to  cycle.  In  many  cases, 
opposite  trends  in  the  spatial  distributions  of  wave 
heights  and  wind  speeds  are  found.  The  statistics  of 
bias,  rms  difference,  linear  regression  coefficients  and 
correlation  coefficient  from  the  32  cycles  of  both 
tracks  are  analyzed  (Figure  8  shows  results  of  3-degree 
averages).  A  rather  large  percentage  (-50%)  of  cases 
show  poor  agreement  based  on  a  combination  of  low 
correlation  and  large  rms  difference  or  bias,  and  the 
regression  coefficient  differs  significantly  from  one. 

The  comparisons  presented  here  indicate  that  our 
numerical  wave  simulation  skills  are  good  in  the 
projection  of  temporal  evolution  but  relatively  poor  for 
the  spatial  variation  of  the  wave  field.  The 
measurements  from  remote  sensing  devices,  including 
spacebome  and  airborne,  are  of  spatial  distribution  in 
nature.  These  data  represent  a  precious  source  for 
enhancing  our  understanding  of  the  spatial  properties 
of  the  wind  and  wave  fields.  The  spacebome  data  are 
especially  useful  for  studying  the  regional  wind  and 
wave  climate,  in  which  case,  the  spatial  density  of  in- 
situ  measurements  is  usually  not  very  high.  Numerous 
comparison  studies  have  shown  that  the  accuracy  of 
altimeter  output  of  significant  wave  height  and  wind 
speed  frequently  exceeds  the  model  output  and  is 
comparable  in  quality  with  in-situ  ocean  buoy 
measurements. 
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Table  1 .  Comparison  ofMean  and  Standard  Deviation  of  WAM  and  TOPEX  Average  Data. 

- 1 - - - tt;  tt  _  _•  r  WinH  fenced 


iaoic 

Latitude/Day 

limit 

Wave  Height 

<//>  (m)  c^J<Hs> 

Wind  Speed 

<L/l0>(m/s)  otn(y'<t^l0> 

TOPEX  1  WAM 

TOPEX  WAM 

TOPEX  WAM  1  TOPEX  WAM 
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0.47 

053 
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0.49 

0,43 

N3  (31.5,  34.5) 

1.51 
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7.07 
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0.50 

0.38 

(b)  Track  69 

D9  (27.5, 36.5) 

1.45 
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0.58 

0.54 

6.66 

6.94 

0.51 

0.49 

S3  (27.5,  30.5) 
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0.46 

0.51 
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1.53 

0.65 
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Track  26  — - 

Q1  (1,90) 
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1.59 
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0.40 
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0.43 

0.36 
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0.63 
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0.36 

8.33 

9.76 

0.29 

0.23 

(b)  Track  69 

01  (1,90) 

1.63 

0.76 

0.48 

0.56 

7.60 

8.10 

0.41 

0.44 

02  (91, 180) 

0.96 

1.07 

0.40 

0.36 

4.50 

5.17 

0.57 

0.38 

Q3  (181,  270) 

1.46 

1.67 

0.69 

0.55 

6.47 

6.74 

0.55 

0.51 

04  (271,360) 

1.69 

1.67 

0.50 

0.50 

7.85 

7.76 

0.38 

0.45 

Table  2.  Statistical  Coefficients  of  WAM  and  TOPEX  Correlations. 


Wave  Height,  //? 

Wind  Speed,  L/m 

Latitude  Limit 

B(m) 

1  A(m)  R 

Cv  c  S(m/s)  A  (m/s)  |  H 

(a)  Track  26 

D9  (25.5,  34.5) 

1.06 

1.07 

0.13 

0.33 

0.93 

1.10 

1.13 

1.07 

1.92 

0.83 

S3  (25.5,28.5) 

1.07 

1.09 

0.17 

0.43 

0.92 

1.02 

1.04 

0.37 

1.73 

0.87 

M3  (28.5,31.5) 

1.07 

1.09 

0.15 

0.41 

0.91 

1.16 

1.21 

1.60 

2.56 

0.78 

N3  (31.5,  34.5) 

1.03 

1.06 

0.06 

0.41 

0.90 

1.09 

1.13 

1.20 

2.52 

0.75 

(b)  Track  69 

D9  (27.5,  36.5) 

1,06 

1.08 

0.08 

0.28 

0.93 

1.06 

1.08 

0.57 

1.33  1 

0.91 

S3  (27.5,  30.5) 

1,05 

1.08 

0.18 

0.42 

0.86 

0.98 

1.01 

0.06 

1.70 

0.87 

M3  (30.5,  33.5) 

0.97 

0.99 

0.04 

0.35 

0.93 

1.00 

1.04 

0.63 

1.97 

0.86 

N3  (33.5, 36.5) 

1.03 

1.10 

0.02 

0.46 

0.83 

1.07 

1.13 

1.01 

2.32 

0.78 

